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Abstract 


In  this  paper,  we  give  an  explicit  method  for  mapping  any  simply  con¬ 
nected  surface  onto  the  sphere  in  a  manner  which  preserves  angles.  This 
technique  relies  on  certain  conformal  mappings  from  differential  geometry. 
Our  method  provides  a  new  way  to  automatically  assign  texture  coordinates 
to  complex  undulating  surfaces.  We  demonstrate  a  finite  element  method 
that  can  be  used  to  apply  our  mapping  technique  to  a  triangulated  geometric 
description  of  a  surface. 


1  Introduction 

The  technique  of  texture  mapping  is  based  on  mapping  an  image  either  synthesized 
or  digitized  onto  a  given  surface.  The  computer  graphics  literature  has  many  such 
works  on  this  topic,  e.g.,  see  [12]  and  the  references  therein.  We  will  not  review 
all  the  literature  on  this  subject  here,  but  just  some  of  the  key  works  in  order  to 
contrast  our  approach. 

Parameterization  of  polygonal  geometric  surfaces  has  direct  use  for  what  Gomes 
et  al.  refer  to  as  image  mapping  [13],  a  category  of  techniques  which  includes  tex¬ 
ture  mapping  [6]  [14]  [15],  reflection  or  environment  mapping  [4],  bump  mapping 
[5],  and  light  mapping  [20].  As  its  name  implies,  image  mapping  associates  a  coor¬ 
dinate  in  image  space  with  every  location  on  the  surface  of  a  geometric  object.  The 
discrete  samples  of  an  image  raster  addressed  by  the  coordinates  in  image  space 
arc  applied  to  the  geometric  object. 

The  mathematical  basis  of  this  mapping  determines  whether  the  image  pro¬ 
vides  the  object  with  the  appearance  of  surface  texture,  a  specular  reflection,  or 
some  other  effect.  In  this  paper,  we  will  focus  only  on  mappings  such  as  texture 
mapping  that  arc  a  function  of  surface  location  and  not  those,  for  example,  based 
on  surface  orientation. 

No  matter  what  kind  of  image  mapping  is  used,  the  parameterization  function 
that  relates  the  geometry  of  the  object  to  image  space  must  be  determined.  Three 
rough  categories  of  mapping  functions  exist  [25]  [13].  The  first  is  a  simple  change 
of  coordinates  from  geometry  space  to  image  space.  Surface  patches  and  other 
parametric  surfaces,  for  instance,  can  often  be  mapped  using  their  natural  paramet¬ 
ric  coordinate  system.  The  second  category,  described  by  Bier  and  Sloan  [3],  uses 
an  auxiliary  surface  as  an  intermediate  between  geometry  and  image.  The  third 
category  projects  the  surface  onto  another  surface  or  geometric  solid.  For  com¬ 
plex  non-parameteric  objects,  finding  a  usable  parameterization  function  can  be  a 
difficult  task. 

In  practice,  the  techniques  used  to  assign  parameters  to  undulating  surfaces  are 
often  ad  hoc  and  may  suffer  from  distortions,  singularities,  or  other  shortcomings. 
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Several  more  rigorous  techniques  have  been  proposed  to  parameterize  complex  sur¬ 
faces.  Bennis  et  al.  describe  a  flattening  and  cutting  technique  [2],  More  recently, 
Lee  et  al.  have  developed  a  multiresolution  approach  based  on  Loop  subdivision 
[18].  Levy  and  Mallet  use  an  algorithm  for  non-distorted  mappings  for  triangle 
meshes  over  cuts  and  other  discontinuities  [19].  This  mapping  method  combines 
planar  graph  theory  with  penalty  functions  to  preserve  right  angles.  The  use  of  har¬ 
monic  maps  for  mapping  surfaces  with  boundary  conditions  is  suggested  in  [10] 
for  multiresolution  analysis  of  meshes.  This  approach  may  have  the  undesirable 
property  of  triangle  flipping. 

In  this  paper,  we  propose  a  method  based  on  a  key  fact  in  conformal  geom¬ 
etry.  Indeed,  one  can  prove  that  any  surface  without  holes  or  self-intersections 
can  be  mapped  conformally  onto  the  sphere,  and  any  local  portion  thereof  onto  a 
disc.  This  mapping,  called  a  conformal  equivalence,  is  one-to-one,  onto,  and  an¬ 
gle  preserving.  Moreover,  one  can  explicitly  write  down  how  the  intrinsic  metric 
on  the  original  surface  is  transformed  and  thus  areas  and  the  geodesics  as  well. 
Specifically,  the  elements  of  the  first  fundamental  form  [E,  F.  G )  are  transformed 
as  ( pE ,  pF.  pG )  for  some  positive  function  p  depending  on  the  point  of  the  sur¬ 
face.  For  this  reason,  conformal  mappings  are  often  described  as  being  similarities 
in  the  small.  This  means  that  the  mapping  can  be  used  to  map  any  planar  image 
or  structure  on  the  sphere  to  our  given  surface  in  an  angle -preserving  manner,  and 
vice  versa.  Locally  shape  is  preserved  and  distances  and  areas  are  only  changed 
by  a  scaling  factor. 

We  should  note  that  the  mapping  we  construct  is  bijective  (onto  and  one-to- 
one),  and  thus  there  is  no  problem  with  triangles  “flipping”  or  overlapping  when 
mapped  to  R2.  Further,  the  method  does  not  require  cuts  to  be  made  on  the  surface. 
The  mapping  is  seamless  and  thus  avoids  the  problem  of  blending  textures  together 
across  cuts. 

Besides  texture  mapping,  one  may  use  this  technique  to  find  longitude  and  lat¬ 
itude  lines  as  well  as  a  north  and  south  pole  on  a  given  surface  S  by  mapping  the 
corresponding  geometric  features  of  the  sphere  onto  S  in  this  conformal  manner. 
Since  angles  arc  preserved,  the  images  of  the  longitudes  and  latitudes  arc  guar¬ 
anteed  to  be  orthogonal  to  each  other.  Surface  deformation  and  flattening  is  also 
very  useful  for  certain  problems  in  medical  imaging  such  as  functional  magnetic 
resonance;  see  [1,  7,  23,  24]  and  the  references  therein. 

What  makes  our  method  practical  and  easy  to  implement  is  the  key  observa¬ 
tion  that  the  conformal  transformation  may  be  obtained  as  the  solution  of  a  second 
order  elliptic  partial  differential  equation  (PDE)  on  the  surface  to  be  deformed. 
For  triangulated  surfaces,  there  exists  a  powerful,  reliable  finite  element  procedure 
which  can  be  employed  to  numerically  approximate  the  conformal  mapping  func¬ 
tion.  The  heart  of  this  procedure  simply  involves  solving  a  pair  of  systems  of  linear 
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equations. 

We  should  note  that  beyond  image  mapping,  the  ability  to  parameterize  ar¬ 
bitrary  objects  of  a  given  topology  has  a  variety  of  applications.  For  example, 
physical  properties  or  other  variables  can  be  associated  uniquely  with  locations  on 
a  surface. 

The  outline  of  this  paper  is  as  follows.  In  Section  2,  we  sketch  the  analytical 
procedure  to  find  the  conformal  mapping.  The  details  of  this  arc  filled  in  the  Ap¬ 
pendix  of  Section  6.  In  Section  3,  we  describe  how  the  numerical  algorithm  works 
on  a  triangulated  surface.  In  Section  4,  we  demonstrate  our  procedure  on  some 
surfaces  both  synthetic  and  digitized,  and  in  Section  5  we  make  some  conclusions 
about  further  research. 

2  Brief  Sketch  of  Mathematical  Theory 

In  this  section,  we  outline  the  mathematical  justification  of  our  surface  mapping 
procedure.  See  [1]  and  the  Appendix  of  Section  6  for  more  details.  We  start  with 
the  basic  assumption  that  the  given  surface  is  a  topological  sphere.  This  means  the 
surface  has  no  holes,  handles,  or  self-intersections. 

Let  S  denote  a  such  a  surface  in  R3.  Let  8P  denote  the  Dirac  delta  (impulse) 
function  atp  E  S,  A  the  Laplace-Beltrami  operator  on  £\{p},  and  i  the  square  root 
of  —1.  (The  Laplace-Beltrami  operator  is  the  generalization  of  the  usual  Laplacian 
operator  to  a  smooth  surface.  See  [21,  11]  and  Section  6.)  Let  S 2  denote  the  unit 
sphere  in  R3  and  let  C  be  the  complex  plane.  Recall  from  the  Introduction  that 
a  conformal  equivalence  is  a  one-to-one,  onto  mapping  which  preserves  angles. 
We  can  now  state  the  following  result  which  provides  the  analytical  basis  for  our 
texture  mapping  procedure: 

A  conformal  equivalence  z  :  E\{p}  — >  S'2\{north  pole}  =  C  may  be  obtained  by 
solving  the  equation 

=  {k  ~  *1) (1) 

The  derivation  of  this  equation  may  be  found  in  the  Appendix  of  Section  6.  This 
result  means  that  we  can  get  the  conformal  equivalence  by  solving  a  second  order 
partial  differential  equation  on  the  surface.  Fortunately,  on  a  triangulated  surface, 
this  may  be  earned  out  using  a  finite  element  technique  we  will  describe  below. 
Also,  we  should  note  that  the  equivalence  above  is  unique  up  to  scaling  and  trans¬ 
lation. 
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3  Finite  Element  Approximation  of  Conformal  Mapping 

We  have  just  sketched  the  analytical  basis  for  conformally  mapping  a  surface  onto 
a  sphere.  We  want  now  to  describe  a  numerical  procedure  for  carrying  this  out.  We 
now  assume  that  E  is  a  triangulated  surface.  We  continue  to  use  the  notation  of  the 
previous  section,  and  let  <7  =  ABC  be  a  triangle  on  whose  face  the  point  p  lies. 

3.1  Finite  Element  Formulation 

We  will  now  briefly  outline  the  finite  element  method  for  finding  the  conformal 
mapping  z  which  is  the  solution  to  the  PDE  (1).  See  [16]  for  details  about  this 
method. 

First  we  construct  a  weak  formulation  of  the  PDE  (1),  i.e.  we  multiply  the 
equation  with  an  arbitrary  smooth  function  /  and  integrate  by  parts  to  get 

/  Lv*-*fis  =  f  l  a -*1)^ 

where  Vz  is  the  gradient  with  respect  to  the  induced  metric  on  E.  One  can  show 
that  z  satisfies  (1)  if  and  only  if  (2)  holds  for  all  smooth  functions  /. 

Next,  we  restrict  our  attention  to  the  finite  dimensional  vector  space  PL(E)  of 
piecewise  linear  functions  on  E,  and  seek  a  z  E  PP(E)  such  that  (2)  holds  for  all 
/  E  PP(E).  For  each  vertex  P  E  E,  let  fp  be  the  continuous  function  such  that 

'  MP)  =  1, 

<  4>p{Q)  =0,  Q  f  P,  Q  a  vertex,  (3) 

(j pp  is  lineal-  on  each  triangle. 

Then  these  fj>  form  a  basis  for  PP(E),  and  so  any  z  E  PP(E)  can  be  written 
as 

z  =  X]  Zp(t)p 

P  vertex  of  E 

for  some  complex  numbers  zp.  Further,  since  (2)  is  linear  in  /,  it  is  enough  to  show 
that  (2)  holds  whenever  f  =  4>q  for  some  Q. 

In  short,  we  want  to  find  a  vector  of  complex  numbers  z  =  (zp),  containing 
one  element  per  vertex,  such  that  for  all  Q, 

J2Zp  J  J  =  ^(p)  -  (4) 
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3.2  Formulation  in  Matrix  Terms 


The  formulation  (4)  is  simply  a  system  of  linear  equations  in  the  complex  un¬ 
knowns  zp. 

Let  us  introduce  the  real  matrix  D  =  ( Dpq  )  where 

Dpq  =  J  j  V(pp  •  S/<pQdS, 

for  each  pair  of  vertices  P,  Q.  It  is  easy  to  see  that  if  P  and  Q  arc  not  connected  by 
an  edge  of  a  triangle,  then  Dpq  =  0.  Thus  the  matrix  D  is  sparse. 

Suppose  that  P  /  Q,  and  PQ  is  an  edge  belonging  to  two  triangles,  say  PQR, 
and  PQS.  A  classical  formula  from  finite-element  theory  [16],  then  says  that 

Ppq  = -^  {cot  ZP  +  cot  ZP}  ,  (5) 

where  ZP  is  the  angle  at  the  vertex  P  in  the  triangle  PQR,  and  IS  is  the  angle  at 
the  vertex  S  in  the  triangle  PQS. 

If  P  =  Q,  then  one  has 

Dpp  =  -  ^2  DPQ ■  (6) 

p^Q 

We  now  consider  the  right-hand  side  of  (4).  Recall  that  p  is  a  point  on  the  face 
of  some  triangle  a  =  ABC.  If  Q  ^  {A,  P,  C}.  then  <pQ  is  identically  zero  on  cr 
and  so  the  right-hand  side  of  (4)  is  zero. 

So  suppose  Q  E  {A,  P,  C},  and  choose  the  u  and  the  v  axis  so  that  A  and  P 
are  along  this  axis,  and  the  positive  v  axis  points  towards  C.  Since  (pq  is  linear  on 
cr, 

dcPg  cPq(B)  -  cPq(A)  dcpQ  <Pq(C)  -  cPg(E) 
du  ||P -A||  ’  dv  \\C-E\\ 

where  E  is  the  orthogonal  projection  of  C  on  AP. 

We  derive  explicit  formulas  for  E  and  rpq(E)  as  follows.  Let  9  be  such  that 

E  =  A  +  9(B  -  A). 

Since  (C  —  E)  T  (P  —  A),  we  have 

{C  -  A-9{B  -  A),B  -  A)  =0, 

so 

(C1  —  A,  P  —  A) 

IIP  -  AN2  ' 


6 


Therefore,  again  by  the  linearity  of  (f>q  on  the  triangle  a,  we  have 

(f>Q(E)^cf>Q(A)+e(cf>Q(B)-(f>Q(A)). 

If  we  define  two  real  vectors  a  =  (oq)  =  (^p(p))  and  b  =  ( bQ )  =  (^-(p)), 
then  the  above  formulas  yield 


aQ  ~  ibQ  :=  < 


0 


-l 


P^W 
P^ 
lW=E\ 1 


+  i 
+  i 


i-e 

\c=m 


Qi{A7B7C}7 
Q  =  A , 

q  =  b7 

Q  =  C. 


and  the  system  (4)  becomes,  in  matrix  terms, 

Dz  =  a  —  ib. 


(7) 


3.3  Solving  the  System  of  Equations 

The  system  of  linear  equations  (7)  can  be  solved  using  standard  techniques  from 
numerical  linear  algebra.  In  practice,  we  solve  for  the  real  and  imaginary  parts  of 
z  =  x  +  iy  separately,  i.e.  we  solve  Dx  =  a  and  Dy  =  b. 

Note  that  since  J2q^pq  =  0  for  all  P,  the  matrix  D  =  (D[>q)  is  singular, 
and  so  we  need  to  show  that  solutions  to  Dx  =  a  and  Dy  =  b  exist.  In  addition, 
we  will  show  that  D  enjoys  several  properties  which  make  these  solutions  easy  to 
compute  numerically. 

We  remark  that  if  Dx  —  0  for  some  non-zero  vector  x  =  (xp).  then  all  the 
elements  of  x  arc  the  same.  To  demonstrate  this,  suppose  Dx  =  0.  Then  clearly 

T,  DpQXpXQ  =  0.  (8) 

P,Q 

Further,  by  definition  of  the  matrix  Dpq  we  have 

[  f  \Vu\2dS  =  DpqXpXq 
7  72  Pq 

where  u  E  PL(S)  is  the  function  with  u(Q)  =  xq  for  all  vertices  Q. 

(8)  and  (9)  together  imply  that  u  is  constant,  and  hence  that  all  xq  are 
conclude  that  the  kernel  of  D  is 

H  :=  {A(1,1,...,1)t  |  A  E  R}. 

By  construction,  D  is  a  sparse,  real,  symmetric  matrix.  The  relation  (9)  implies 
that  D  is  positive  semi-definite,  and  together  with  the  analysis  above,  we  see  that  D 


(9) 

Equations 
equal.  We 
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maps  if-1,  the  orthogonal  complement  of  H,  bijectively  to  itself.  Thus  the  equation 
Dx  =  f  is  solvable  if  and  only  if  /  E  HL  (i.e.  Jq  —  0)’  and  this  solution  is 
unique  up  to  addition  of  an  element  of  H.  We  note  that  a  and  b  are  indeed  in  H±. 

Since  D,  restricted  to  II  is  sparse,  symmetric  and  positive  definite,  equations 
Dx  =  a  and  Dy  =  b  are  particularly  well  suited  for  numerical  solution  by  the  well- 
known  conjugate  gradient  method.  Although  D  is  singular,  this  method  involves 
only  multiplications  by  D  and  addition  of  vectors  in  H2-,  and  so  quite  literally 
solves  the  equations  for  D  restricted  to  II ~ . 

The  speed  of  the  conjugate  gradient  method  allows  us  to  find  x  and  y  in  a 
combined  time  of  under  3  seconds  for  a  surface  with  8128  triangles,  using  a  Sun 
UltraSparc  10. 

3.4  Algorithm  for  Texture  Mapping 

We  may  summarize  our  procedure  for  texture  mapping  via  the  construction  of  a 
conformal  map  z  as  follows: 

Input:  A  triangulated  surface  E,  and  an  image  to  be  used  as  a  texture. 

(1)  Compute  the  elements  of  the  matrix  D ,  and  the  vectors  a  and  b,  using  the 
formulas  from  section  3.2. 

(2)  Solve  the  systems  of  linear  equations  Dx  —  a  and  Dy  =  —b. 

(3)  Compose  z  =  x  +  iy  with  inverse  stereo  projection  to  get  a  conformal  map 
to  the  unit  sphere.  Specifically,  send  the  point  x  +  iy  to  the  point  {2x/{l  + 
r2),2y/(l  +  r2),2r2/(l  +  r2)  —  1),  where  r2  =  x2  +  y2. 

(4)  Compute  the  spherical  coordinates  of  the  points  on  the  sphere.  These  will  be 
new  texture  coordinates. 

Output:  Texture  coordinates,  and  rendered  original  surface,  using  the  new  texture 
coordinates  to  index  into  the  input  image. 

4  Computer  Experiments 

We  illustrate  the  texture  mapping  procedure  on  a  number  of  surfaces  both  synthetic 
and  digitized. 

Figures  1  through  3  arc  concerned  with  mapping  textures  onto  synthetic  sur¬ 
faces  which  contain  highly  concave  and  convex  areas.  In  Figure  1,  we  conformally 
map  a  mandrill  image  onto  a  synthetic  surface.  Similarly  in  Figure  2,  we  con¬ 
formally  map  a  thermal  image  of  the  earth  onto  a  synthetic  tooth-shaped  surface. 


Notice  how  in  both  cases  the  local  geometry  is  very  well-preserved,  and  the  map¬ 
pings  arc  performed  without  any  tears  or  seams  in  the  image.  Each  image  in  Figure 
3  shows  two  copies  of  the  earth,  one  mapped  to  the  outside  and  another  to  the  in¬ 
side  of  a  vase-shaped  object.  The  rim  of  the  vase  is  not  a  boundary,  but  rather  a 
fairly  sharp  edge  between  the  inside  surface  of  the  vase  and  its  outside.  The  vase 
is  thus  a  topological  sphere.  This  example  demonstrates  a  more  general  property 
of  our  method;  we  have  found  it  to  be  quite  robust  in  the  presence  of  edges  such  as 
those  found  on  polyhedral  surfaces. 

In  Figure  4,  we  performed  a  white  matter  segmentation  of  the  brain  whose 
points  arc  colored  according  to  mean  curvature.  (For  details  about  how  such  a 
segmentation  may  be  performed  we  refer  the  reader  to  [17,  23].)  To  get  a  surface 
triangulation,  we  used  the  Visualization  Tool  Kit  (VTK).  Using  our  conformal  geo¬ 
metric  technique,  we  then  map  the  highly  convoluted  brain  surface  onto  the  sphere. 
Notice  how  the  sulci  (indentations)  and  gyri  (protrusions)  are  clearly  represented 
on  the  sphere  with  a  nice  preservation  of  the  local  geometry.  Such  mapping  tech¬ 
niques  have  uses  in  functional  magnetic  resonance  imaging  for  the  brain  flattening 
problem.  Since  it  is  important  to  visualize  functional  MRI  data  for  neural  activity 
within  the  folds  of  the  brain  surface,  flattening  methods  have  become  an  important 
area  of  research  in  medical  imaging;  see  [1,  7,  23,  24]  and  the  references  therein. 

Figure  5  shows  how  orthogonal  coordinate  systems,  i.e.  parameterizations, 
can  be  established  on  surfaces  using  our  technique.  The  colored  areas  correspond 
to  spherical  rectangles  with  boundaries  lying  on  lines  of  longitude  and  latitude. 
Notice  that  when  these  rectangles  are  mapped  from  the  sphere  to  these  surfaces,  the 
corners  of  the  rectangles  continue  to  meet  at  right  angles,  due  to  the  conformality  of 
the  mapping.  As  mentioned  before,  the  conformal  mapping  from  the  surface  to  the 
sphere  is  not  unique.  Essentially,  this  is  because  we  can  compose  the  map  with  any 
conformal  mapping  from  the  sphere  to  itself.  However,  if  we  specify  the  image  on 
the  sphere  of  three  points  on  the  original  surface,  then  there  is  only  one  conformal 
mapping  which  satisfies  this  condition.  We  can  take  advantage  of  this  fact,  because 
it  allows  us  to  specify  which  two  points  on  the  surface  will  correspond  to  the  north 
and  south  poles  on  the  sphere.  In  Figure  5,  we  have  chosen  these  points  so  that  the 
south  pole  is  as  far  away  as  possible  from  the  north. 

Figure  6  shows  how  decimation  of  triangles  affects  the  mapping.  One-fourth 
of  the  triangles  on  each  surface  were  removed  to  create  the  next  surface,  and  the 
conformal  mapping  was  re-calculated.  From  left  to  right  the  surfaces  arc  composed 
of  32512,  8128,  2032  and  508  triangles.  We  made  consistent  choices  for  north 
and  south  poles  on  each  surface  so  that  the  images  could  be  compared  visually. 
Removal  of  triangles  tends  to  flatten  out  portions  of  the  surface,  and  this  change  in 
geometry  results  in  a  slight  shifting  of  the  texture. 

Figure  7  shows  the  effect  of  affine  transformations  on  the  procedure.  The  sur- 
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face  in  the  upper  left  comer  was  transformed  three  times  by  multiplying  its  vertices 
by  fixed  matrices  chosen  at  random.  Conformal  mappings  to  the  sphere  were  then 
calculated.  Again,  the  poles  were  chosen  consistently  to  aid  in  visual  comparison. 
Affine  transformations  alter  surface  angles,  and  so  naturally  the  texture  coordinates 
we  calculate  arc  not  completely  affine  invariant.  However,  we  note  that  continu¬ 
ously  deforming  the  surface  results  in  continuous  changes  in  texture  coordinates 
and  so  singularities  in  the  texture  map  do  not  form.  Here,  although  the  differences 
between  the  surfaces  arc  rather  severe,  the  texture  mapping  remains  relatively  con¬ 
sistent.  Given  these  results,  we  arc  very  encouraged  that  this  process  may  be  used 
for  automatic  texture  mapping.  A  simple  example  would  be  applying  a  mottled 
apple-skin  texture  to  any  surface  that  is  reasonably  apple-shaped. 

5  Conclusions 

In  this  section,  we  summarize  a  number  of  the  key  aspects  of  our  method.  Fur¬ 
ther,  we  show  how  the  method  avoids  some  of  the  problems  often  associated  with 
surface  mapping,  and  discuss  applictions  and  areas  for  possible  future  research. 

5.1  General  Remarks 

In  this  paper,  we  have  described  a  general  method  for  mapping  certain  surfaces  to 
the  sphere  in  a  manner  which  preserves  the  local  geometry.  Since  angles  are  pre¬ 
served  by  this  mapping,  a  texture,  when  applied  to  the  surface,  has  much  the  same 
appearance  that  is  has  in  the  plane  or  on  the  sphere.  The  method  also  provides  an 
effective  way  to  establish  orthogonal  coordinate  systems  on  these  surfaces.  Since 
the  mapping  is  naturally  bijective,  there  is  no  inherent  problem  with  triangles  over¬ 
lapping  or  “flipping”  when  mapped  to  the  plane.  Nor  does  the  method  require  cuts 
to  be  made  on  the  surface;  it  is  a  seamless  mapping  and  so  avoids  the  problem  of 
blending  together  textures  across  cuts.  The  formulation  of  the  partial  differential 
equation  which  is  central  to  the  method  provides  it  with  an  attractive  simplicity  and 
elegance. 

The  simplicity  of  the  equation  formulated  allows  us  to  use  standard,  powerful 
numerical  methods  to  solve  it.  The  finite  element  method  from  numerical  partial 
differential  theory  applied  to  our  system  with  the  proper  boundary  conditions  for¬ 
mulated  above,  together  with  the  conjugate  gradient  method  from  numerical  linear 
algebra,  form  a  robust  and  reliable  solution  system.  As  seen  in  the  examples,  this 
robustness  shows  up  in  the  stability  of  the  method  when  applied  to  surfaces  of 
varying  shape  and  triangulation,  even  in  the  presence  of  sharp  edges.  In  addition, 
the  speed  of  the  conjugate  gradient  method  makes  the  finding  of  the  solution  quite 


10 


practical. 


5.2  Other  Applications 

We  have  shown  how  our  method  can  be  useful  in  automatic  texture  mapping,  and 
in  the  construction  of  coordinate  systems  on  even  highly  undulated  surfaces.  We 
have  already  discussed  the  application  of  our  technique  in  medical  imagery  for 
functional  magnetic  resonance  data  in  order  to  solve  in  a  natural  manner  the  brain 
flattening  problem.  Moreover,  a  number  of  pathologies  have  been  associated  with 
deformations  of  brain  structures.  We  are  very  hopeful  that  our  conformal  method¬ 
ology  will  be  useful  in  quantitatively  describing  such  pathologies. 

Our  conformal  geometric  approach  can  also  be  utilized  for  several  other  prob¬ 
lems  in  computer  graphics  including  defining  the  north  and  south  poles  on  a  given 
surface  and  in  the  computation  of  shading  maps.  Traditional  computer  graphics 
shading  is  done  by  applying  synthetic  lighting  calculations  at  each  polygon  vertex, 
then  interpolating  between  the  vertices  to  find  approximations  to  those  values  on 
the  interior  of  the  polygons.  The  cost  of  this  calculation  obviously  goes  up  with  the 
number  of  polygons.  However,  coarsely  tesselated  objects  appeal-  more  realistic  if 
shading  is  done  at  a  finer  scale,  and  finely  tesselated  objects  are  expensive  to  render 
in  part  because  of  the  large  number  of  lighting  calculations.  Also,  it  is  harder  to 
do  complex  shading  operations  in  hardware  where  vertex-level  shading  is  usually 
performed.  Shading  maps  get  around  several  of  these  problems  in  a  number  of 
applications.  Given  a  mapping  from  3D  geometry  to  a  2D  texture,  lighting  and 
interpolation  can  be  performed  on  a  regular  grid  in  two  dimensional  space  instead 
of  three  dimensions.  Texture  mapping  is  then  used  to  map  the  result  of  the  shad¬ 
ing  calculations  back  onto  the  surface.  We  plan  to  apply  our  methodology  to  this 
problem  as  well. 

5.3  Future  Research 

We  are  currently  working  on  the  problem  of  global  area  and  length  distortion.  It 
is  a  classical  result  in  differential  geometry  (see,  for  example,  [8]),  that  it  is  not 
possible  to  map  a  general  surface  to  the  sphere  or  plane  without  some  distortion  of 
lengths,  and  so  at  best  a  compromise  must  be  made  between  length  preservation 
and  conformality  of  the  map. 

We  are  testing  certain  variational  methods  to  find  a  suitable  compromise  be¬ 
tween  our  conformal  mappings  and  ones  with  minimal  area  and  length  distortions. 
We  are  also  working  on  extending  our  method  to  surfaces  with  more  general  topol¬ 
ogy,  such  as  surfaces  with  holes  or  boundaries. 
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6  Appendix:  Derivation  of  Partial  Differential  Equation 


In  this  section,  we  outline  the  derivation  the  partial  differential  equation  (1).  See 
also  [1,  11,  21]  and  the  references  therein.  Let  p  be  a  point  on  the  surface  £, 
z  :  £  — >  S'2  be  a  conformal  equivalence  which  sends  p  to  the  north  pole,  and 
let  ( u ,  v)  be  conformal  coordinates  on  £  near  p  such  that  u  =  v  =  0  at  p. 
Conformal  coordinates  u.  v  are  such  that  the  metric  at  the  point  p  is  of  the  form 
ds 2  =  X (u,  v)2(du 2  +  dv2).  We  can  always  ensure  that  at  the  particular  point  p, 
X (p)  =  1,  VX(p)  =  0.  It  can  be  proven  that  such  conformal  coordinates  always 
exist  [8],  In  these  coordinates,  the  Laplace-Beltrami  operator  takes  the  simpler, 
more  familial-  form 


1  /  d 2  d2  \ 

X(u,v)2  yriu2  dv2  J  ' 

Set  w  —  u  +  iv.  Since  the  mapping  z  =  z(w)  is  one-to-one,  it  follows  that  z 
has  a  simple  pole  at  w  =  0,  and  hence  a  Laurent  series  representation 

A 

z(w)  = - 1-  B  +  Cw  +  Dw 2  +  . . . 

w 

We  interpret  this  equation  in  the  sense  of  distributions  and  apply  A  to  both  sides  to 
get 

AZ  =  (w)  ’ 

the  non-negative  powers  of  w  being  harmonic  near  p.  We  need  only  find  z  up  to  a 
constant  multiple,  so  taking  A  = 

^a(£)  =  (l;_il;)losW 

1  f  d  .  d  N 

=  iU-ri"051”1 

1  /  d  d  \ 

where  we  have  used  the  fact  that  ^  log  |u’|  is  the  fundamental  solution  for  the 
operator  A.  This  gives  (1), 
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Figure  1 :  Mandrill  Image  on  Synthetic  Surface 
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Figure  2:  Thermal  Image  of  Earth  on  Synthetic  Surface 
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Figure  3:  Earth  on  Inside  and  Outside  of  a  Vase 


Figure  4:  Two  Views  of  Flattened  White  Matter 


Figure  5:  Orthogonal  Coordinate  Systems 


Figure  7 :  Effect  of  Affine  Transformations 
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